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ABSTRACT 

The geometry of the dust distribution within the inner regions of Active Galac- 
C"| . tic Nuclei (AGN) is still a debated issue and relates directly with the AGN unified 

O [ scheme. Traditionally, models discussed in the literature assume one of two distinct 

dust distributions in what is believed to be a toroidal region around the Supermas- 
sive Black Holes: a continuous distribution, customarily referred to as smooth, and a 
concentration of dust in clumps or clouds, referred to as clumpy. 

In this paper we perform a thorough comparison betw een two of t h e mos t popular 
models in the lite rature, namely the sm ooth models by iFritz et all (|2006( l and the 
clumpy models bv lNenkova et al.l ()2008aD . in their common parameters space. Partic- 
ular attention is paid to the silicate features at ~9.7 and ~18 pm, the width of the 
infrared bump, the near-infrared index and the luminosity at 12.3 pm, all previously 
reported as possible diagnostic tools to distinguish between the two dust distributions. 
We find that, due to the different dust chemical compositions used in the two models, 
the behaviour of the silicate features at 9.7 and 18 pm is quite distinct between the 
two models. The width of the infrared bump and the peak of the infrared emission 
can take comparable values, their distributions do, however, vary. The near-infrared 
index is also quite different, due partly to the primary sources adopted by the two 
models. Models with matched parameters do not produce similar SEDs and virtually 
no random parameter combinations can result in seemingly identical SEDs. 
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1 INTRODUCTION 

Active galactic nuclei (AGN) can be classified in a variety 
of ways depending on the characteristics of their spectral 
energy distribution (SED) at various wavelengths. The com- 
monly used division in type 1 and type 2 AGN is primarily 
based on their properties in the UV/optical wavelengths, 
with type 1 objects typically showing broad emission lines, 
while type 2 only having narrow emission lines in their spec- 
tra. According to the unified scheme for AGN (|Antonuccil 
1 19931 : lUrrv fc Padovanilll995l h the differences between the 
type 1 and type 2 AGN are an orientation effect, as first 
suggested bv lRees et al.l (| 19691 ). For certain lines of sight the 
dust can obscure the central engine giving rise to a number of 
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differences in the observed SED of AGN at almost all wave- 
lengths. This happens as dust grains, present in an optically 
thick region surrounding the nucleus, absorb ultraviolet pho- 
tons coming from the central region and re-radiate them in 
the infrared (IR). Therefore, the presence of dust around 
the central region of AGN is the key to understanding the 
differences between type 1 and type 2 objects. 

The variety of radiative transfer models developed to 
reproduce the observed dust emission in the IR can be 
divided in two classes: "smooth" models characterised by 
a continuous dust distribution in a toroidal or fiared-disk 
shape in which the density can only vary smoothly within 
the torus, and "clumpy", in which the dust is distributed 
in clumps or clouds. Clumpy models are a more likely rep- 
resentation of the real dust distribution as a smooth dust 
distribution would result in collisions that would raise the 
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tem perature to levels too hig h for the dust to survive (see 
e.g. iKrolik fc Begelmanlll988r ). On the other hand, smooth 
models were the first to be developed, being computation- 
ally simpler and, in many aspects, a good approximation 
when calculating the 1R SED of AGN. 

Smooth models, however, confront important issues 
when used to match the observations: the IR bump 
they produced was narrower t han what was observed 
IIDullemond fc van Bemmell [20051 . see also ij3.2[) ; the 10 um 
silicate feature f i)3.1|) was often observed in absorption in 
type 2 sources, but had ba rely been seen in emission in either 
type 1 or type 2 sources jDullemond fc van Bemmell 120051 ) 
until recent observations with Spitzer/IRS; and the same 
feature in absorption in type 2 views had always been ob- 
served to be shallower than what t he model predicted. With 
this in mind, iNenkova et"afl (|2002h were the first to present 
a torus model with dust distributed in clumps, with the sil- 
icate emission feature attenuated and the IR bump broad 
enough to fit the observations. The same models, however, 
fail to reproduce the sh ort wavelength e mission emerging 
from the hot dust (e.g. iMor et all 120091 ) in type 1 AGN. 
This behaviour of the clumpy models is not specific to this 
particular set of models but appears in other clumpy model 
realisations fe.g. IPolletta et alll2008l ). 

The success of both classes of models in fitting dif- 
ferent parts of the observed AGN SEDs keeps the issue 
of the dust distribution in AGN open, as no conclusions 
can be drawn from the simple comparison between observed 
and model SEDs. The differences arising from the two ge- 
om etries have so far been the focus of two different stud- 
ies. iDullemond fc van Bemmell |2005l ) compared their own 
two-dimensional radiative transfer models of smooth and 
clumpy tori in terms of the resulting width of the SED, 
strength of the silicate feature at 9.7 um and isotropy pa- 
rameter. They concluded that, despite the distinct nature 
of the models and the variations this may cause to the 
shape of the SEDs, distinguishing between the two distri- 
butions based on t he br oad band SEDs was not possible. 
ISchartmann et all |2008T l implemented a three dimensional 
clumpy mo del and compared it wit h their previous continu- 
ous model l|Schartmann et al"1l2005l ) as well as other clumpy 
configurations. Their analysis, with an emphasis on the be- 
haviour of the 9.7 um silicate feature, confirmed their pre- 
vious conclusion that the mid-IR SEDs of AGN are mainly 
determ ined by the innermos t part of the torus. More re- 
cently, IStalevski et al.l (|2012l ) explored the implications of 
clumpy and continuous dust distributions on the IR SEDs, 
by expl oiting the three-d imensional radiative transfer code 
SKIRT |Baes et al.ll201ll ). and found globally no significant 
dissimilarities in their full set of models. 

Other than the works mentioned above, a systematic 
comparison between the models currently used in the litera- 
ture is missing. The work we present here intends to partly 
fill this gap by providing a thorough comparison between the 
two probably most used smooth and clumpy models in the 
literature to date, namely an updated grid of smooth models 
based on Fritz, Franceschini & Hatziminaoglou (2006) (see 
32. II f or the de t ails) an d the updated clumpy model grid by 
INenkova et al.l l|2008al ). The two sets of models and the pa- 
rameter space they cover are described very briefly before 
constructing restricted grids of the two classes of SEDs by 
matching the model parameters (©. We then characterize 



and compare the two sets of SEDs and investigate the de- 
rived distributions of the respective model parameters (Sj3j). 
We finally summarise our results and expose the implica- 
tions for the use of the two approaches of radiative transfer 
models to reproduce the observed AGN SEDs (§3). 



2 DESCRIPTION OF MODELS AND GRIDS 
SELECTION 

In this section we briefly summarize the main features of the 
models compared in this work, describing their main features 
and parameters, as well as the selection of the final grids of 
models on which we perform the compa rison . For the mod- 
els details we refer to iFritz et al l ([20061 ) and INenkova et al.l 
(|2008al lbh. hereafter F06 and N08, respectively. What fol- 
lows is a brief description of the two different model grids. 
For a complete and detailed description see the respective 
papers. 

2.1 Smooth dust distribution 

The grid of smooth models used for this study is based on the 
F06 models, with some minor updates, as described below. 
The original models are among the most popular smooth 
models in the literature to date (see e.g. Rodighiero et all 

' '20091. 



dig hicrc 

|2007l ; iBerta et al I [20071; iHatziminaoglou et all I2OO8I 
| 2010l;lAgol et alJl2009l:|Pozzi et al.ll2010l ; lNatale et all 



2010; 



Vignaliet all 120091 , l201ll ). The mixture of graphite (53%) 



and silicate grains (47%), the distribution of grain sizes 
(Mathis, Rumpl & Nordsieck 1 997) and absorption a nd scat- 
tering coefficients (taken from lLaor fc Drainelll993l ) remain 
unchanged with respect to F06. The spectral index for the 
power la ws describing the central source has been updated, 
following ISchartmann et al. I (|2005l ): 
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A new model grid with a finer sampling of the parameter 
space and a better wavelength resolution has been created. 
Given that the dust density is described in polar coordinates 
as 



p{r,9) — pa ■ r q ■ e 



-7X |cos(0)| 



(2) 



the following parameters and their respective values are ex- 
plored: 

- the torus amplitude, defined as the angular region occu- 
pied by the torus dust, complementary to the opening angle 
of the torus, 6: 60°, 100° and 140°; 

- the parameters of the dust distribution, namely q : 0.00, 
0.25, 0.50, 0.75 and 1.0; and 7: 0.0, 2.0, 4.0 and 6.0; 

- the equatorial optical depth at 9.7 um, r e9 (9.7): 0.1, 0.3, 
0.6, 1.0, 2.0, 3.0, 6.0 and 10.0; 

- the outer-to-inner radius ratio, Y: 10, 30, 60, 100 and 
150. 

The model SEDs are computed at different lines of sight 
with respect to the torus equatorial plane in order to account 
for both type 1 and type 2 object emission, from 0° to 90° 
in steps of 10°. 
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2.2 Clumpy dust distribution 



The N08 cl umpy models, extensively used in the litera- 
ture ( see e.g.lMor et al.ll2009l ; iNikutta et al.ll2009l ; lDeo et all 
l201ll ; iRamos-Almeida et all [20111 ; lEsquei et al.ll2012l ), as- 
sume a grain composition of 53% silicate and 47% g raphite, 
with t he optical constants for the gr aphite taken fromlDrainej 
l|2003h and that of the silicates from lOssenkopf et al, l l| 19921 ). 
The primary source is described by a p i ecewi se power-law 
distribution following I Rowan- Robinson! (| 19951 ), that when 
expressed in terms of L(A) it takes the form: 
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The angular distribution of clumps is a Gaussian of width 
a, given by: 



N T {P) = N e 



(-r/O 



(4) 



where /3 (=90-0 using the F06 notation) is the angle with 
respect to the torus axis, iVo is the average number of clouds 
along radial equatorial rays, where the clouds follow a Pois- 
son distribution. Model parameters and their respective val- 
ues can be summarised as follows: 

- the width of a Gaussian angular distribution of the 
clouds, a, ranging from 15° to 70°, determining the spatial 
distribution; 

- the outer-to-inner radius ratio of the cloud distribution, 
Y= 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 150, and 200; 

- the average number of clouds along a given radial direc- 
tion within the equatorial plane, iVo, taking values between 
1 and 15; 

- the power-law index determining the radial distribution 
of clouds, q between 0.0 and 3.0; 

- the optical depth of a single cloud tv =5.0, 10.0, 20.0, 
30.0, 40.0, 60.0, 80.0, 100.0, 150.0, calculated in the V band 
at 0.55 um. 

Again, the SEDs are created for different viewing direc- 
tions with respect to the equatorial plane, from /3 = 0° to 
P = 90°, in steps of 10°. 



2.3 Intrinsic differences between the two dust 
distributions 

The two types of models, though both able to reproduce a 
variety of observations, are intrinsically quite different. Fig. 
[1] shows an example of a smooth (left) and a clumpy (right) 
model SEDs with matched model parameters (see Sec. 12. 4[) . 
viewed at different angles ranging from 0° to 90° degrees 
from the equatorial plane. The sudden jump of the smooth 
SEDs occurs at the angle where the dust starts intercepting 
the line of sight (70° in this particular case). The clumpy 
SEDs present a smooth transition from one viewing angle 
to the next, due to their Gaussian angular distribution. 

As already mentioned, while the F06 models use the 
silicates absorption and scattering coefficients given in 
lLaor fc Drain? (119931), N08 make use of the values given 
in lOssenkopf et al.l ( 19921) . And while the absorption coeffi- 
cient for the former peaks at A ~ 9.5 um, that of the latter 
peaks at 10 um. In order to keep the notation of the paper 
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Figure 1. Example of smooth (left) and clumpy (right) model 
SEDs with comparable model parameter values, viewed at differ- 
ent inclinations. 



simple but also consistent with other works in the literature, 
we will refer to the "9.7 um silicate feature" throughout but 
the reader should keep this inaccuracy in mind. 

The most important implications of the AGN unified 
scheme, is that type 1 objects are observed along a dust free 
line of sight, while in type 2 objects the dust intercepts the 
view to the nucleus. If the dust has a continuous distribution 
the chance of seeing the central source is a mere question of 
the line of sight. In a clumpy medium the chance to have 
a direct view of the central source is in fact the probability 
to encounter, on average, zero clouds, that is less than 1 
even for edge-on lines of sight. To simplify the study, we 
only consider the two extreme inclinations i.e. edge on (9 = 
90°, Eq. [2| and face-on (8 = 0°) for smooth models, while 
we consider as type 1 and 2 clumpy models those with a 
probability greater and lower than 0.5, respectively, to see 
directly the AGN, for each model parameters combination. 



2.4 The matched model parameter spaces 

The two original model grids cover a large parameter space 
and produce SEDs with a considerable overlap when their 
parameters are "matched" . The selection of the matched pa- 
rameter space is of interest here, as is not always straight 
forward given the different nature of the two models. How- 
ever good analogies can be found, under a few reasonable 
assumptions. The inner-to-outer radius ratio, Y , is a param- 
eter in common between the two models. The radial varia- 
tion of the dust density or that of the clouds, q, can also 
be considered as equivalent between the two models. The 
relation between the optical depth at 9.7 um and that in 
the V band is given by N08 as rg.r = 0.042 x tv- The equa- 
torial optical depth T9.7 for clumpy models is then equal 
to 0.042 x tv x No. The resulting Tg.7 are not identical to 
those used in the smooth models, the values, however, are 
very close and the differences they introduce in the resulting 
SEDs negligible. Finally, considering only the two extreme 
lines of sights, the torus opening angle is of no relevance and 
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SMOOTH 


CLUMPY 
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10, 30, 60, 100, 150 


10, 30, 60, 100, 150 


q 


0, 1 


0, 1 
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2, 4, 6 




IT 




60°, 45°, 35° 


Tg.7 


0.3, 0.6, 1, 2, 3, 6, 10 




No 




1 - 15 


T V 




5, 10, 20, 30, 40, 






60, 80 



Table 1. Matched model parameter considered in the compara- 
tive study. 
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Figure 2. The range of SEDs covered by the smooth (left) and 
clumpy (right) dust configurations in the restricted parameter 
grids. The coverage is shown for both type 1 (filled regions) and 
type 2 (dashed) inclinations. 



7 (Eq. [2j and a (Eq. [4} can be taken such that the distribu- 
tion of dust and clumps match each other (see Table 1 for 
the exact values). 

Table [TJsummarizes the values of the models parameters 
that will be considered henceforth. Fig.[2]shows, for illustra- 
tion purposes, the shapes of the model SEDs characterised 
by the above parameters, for type 1 and type 2 views. 

In matching the models parameters as described above 
we significantly restrict the model grids and end up com- 
paring a total of 614 smooth and 480 clumpy models. The 
difference in the numbers is due to the fact that often more 
than one combinations of No and tv correspond to the same 
Tg.7. Furthermore, since we only consider the two extreme 
inclinations (smooth) and probabilities (clumpy), as already 
explained before Q becomes of no relevance and hence three 
values shown in Table [T] are consider for each given 7 and a. 

Note that in matching the models parameters, some 
parameter values were left out from both model grids and 
the comparison presented here only applies to the restricted 
grids. In effect what is left out from these grids are clumpy 
models with very compact configuration (q > 1) since there 
are no available equivalent smooth models, smooth models 
with very low optical depth (rg.7 < 0.3), as well as clumpy 




1 



. No 



Figure 3. The values of S9.7 for matched smooth and clumpy 
models. The points are colour-coded based on the value of Nq. 



models with various combinations of iVo and tv that do not 
correspond to any 79.7 from the smooth grid. 



3 MODEL-TO-MODEL COMPARISON 

To compare the characteristics of the SEDs obtained from 
the two different models of dust distribution we restrict 
ourselves to the common model grid discussed in fc|2.4l and 
measure characteristic quantities, namely the prominence 
or strength of the silicate feature, S, the width of the IR 
bump, Wir, the near-IR slope, am, the peak of the emis- 
sion, Xp ea k, and the monochromatic luminosity at 12.3 jim. 
We normalise all models to a common accretion luminosity 
of 10 45 erg/sec. 



3.1 The silicate features 

The silicate emission feature at ~9.7 um has been under 
the spotlights because its prominent appearance, mostly in 
smooth models, did not match any observations for many 
years, questioning the reliability of the Unified Scheme and 
the accuracy of the models themselves. In clumpy models, 
on the other hand, this feature is often smeared out in axial 
viewing in spite of its prominence in emission from directly 
individual clouds. The controversy on the prominence of the 
feature was part ly solved by observations carried out with 
Spitzer /IRS (e g. isiebenmorgen et al.|[2005l ; lHao et alJlioOrj ; 
IShi et alj|200a ). 

The strength of the silicate feature is defined as the log- 
arithm of the ratio between the flux F measured within the 
line profile over the continuum flux F c at such wavelength, 



ln(F(A m )/F c (A m )) 



(5) 



where A m is the wavelength at which the feature's strength 
is an extremum with a value in the inte rval between 8.5 an d 
11.5 um (for the computation of F c see lSirockv et al] [2008). 
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As already mentioned in Sec. 12,31 F06 and N08 models con- 
sider different absorption coefficients that peak at different 
wavelengths, but for simplicity we call the strength of the 
silicate feature around 9.7 um S9.7, irrespective of A m . 

Fig. compares the values of Sg.7 for matched smooth 
and clumpy models in a face-on inclination. The symbols 
are colour-coded according to the value of No. For small 
values of No the values of the feature for the two models 
lie close to the 1:1 line (shown in black), however as No 
takes larger values the points deviate from the line, with 
the clumpy models showing a weaker feature. This is due 
to increasing attenuation with increasing number of clouds, 
as dust is optically thick to itself. The distribution of S9.7 
for the two dust configurations can be seen in the two top 
histograms in Fig.[4l with the green (red) lines corresponding 
to smooth (clumpy) configurations. 

Both Fig. [3]and the histograms at the top of Fig.[4]show 
that a large interval of 59.7 values is covered by both dust 
configurations. The clumpy models with type 1 views can 
produce equally strong silicate features in emission (Fig. U 
top left histogram), despite repeated claims in the literature 
to the contrary. A large fraction of them, however, does ex- 
tend to weaker Sg.7 values. The striking differences in the 
behaviour of Sg.7 (with means of 0.79 and 0.62 for type 1 
smooth and clumpy, respectively, and -0.48 and 0.06 for type 
2 views) have their origin in both the different chemical com- 
positions used by the two models (see Sec. 12. 3p and the fact 
that only the restricted model grids are being compared, 
leaving out a number of models (both smooth and clumpy) 
with parameters that can not be matched by the other dust 
distribution. On the other hand, the long standing issue of 
very deep silicate absorption produced by smooth models is 
confirmed, with the tail of the S9.7 distribution for edge-on 
inclinations (Fig. 01 top right histogram) extending to large 
negative values. Both models can produce silicate emission 
in type 2 views. Since objects with such chara cteristics are 
rather uncommon (but not unheard of, see e . g. | Sturm et alJ 
20061; iTeplitz et a"fll2006l ; iMason et all l2009t iNTkutta et all 



2009), the models that produce such features should also be 



seen as non-representative, yet realistic. Furthermore, sili- 
cate in absorption in type 1 views is also produced by both 
models, although only marginally. 

Recently the attention has turned towards the 18 um 
silicate feature which, being broader and fainter, had almost 
escaped detection. Its strength, Sis, is also computed from 
Eq. [5] with its extremum taking values in the interval be- 
tween 17.0 and 19.5 um. Sg.7 and Sis are shown in Fig. [3] 
for type 1 (left) and type 2 (right) views. 

Clumpy models show stronger Sis emission in type 1 
views and span a larger range compared to the smooth 
models. For type 2 views, the values for Sis emerging from 
the two dust configurations are somewhat more similar, but 
again clumpy models extend to stronger features in emission 
while there are a few smooth models that show very deep 
abs orption. 

iThompson et all |2009l ) proposed that the comparison 
of the strengths of the two silicate features could be used as 
a diagnostic of dust composition and, also, to discriminate 
between smooth and clumpy distributions. Fig. [4] shows in- 
deed that the relative strength of the two features is very 
distinct between the two models, however as pointed out by 
ISirockv et al.l (|2008h this difference reflects the effect of the 



different chemical compositions used in t he two models (dis- 
cussed in Sec. 12. 31 see also Figs. 7 and 9 in lSirockv et al 120081 
showing the changing slope in the Sis vs S9.7 distribution 
as a function of the dust chemistry). 



3.2 Infrared SED 

We now examine the properties of the IR model SED result- 
ing from two different dust distributions, namely: 

- the width of the IR bump, Wir, defined as the logio 
of the frequency range in which the spectrum is more 
than 1/3 of its peak va lue (expressed in F v ) as in e.g. 
iGranato fc Danesd (|l994T ); 

- the wavelength where the flux value (in F v ) of the SED 
peaks, \ P eak\ 

- the spectral index at near-IR wavelengths, defined as 



am 



log 10 (F 4 .s) -log 10 (F 3 . 5 ) 



logio(A4.s) -log 10 (A 3 .5) 
the monochromatic luminosity at 12.3 um, L12.3. 



(6) 



The distribution of the first three parameters for type 1 and 
2 views are shown in Fig. [5] 

While the mean of the distributions of Wir is the 
same for smooth and clumpy models (12.88 for type 1 and 
12.9 type 2 views), and the medians differing only slightly 
(12.86 versus 12.88 for type 1 views and 12.89 versus 12.9 
for type 2 views), clumpy models produce, on average, 
wider IR bumps, as already noted by other author s (e.g. 
iDullemond fc van Bemmelll2005l : iNenkova et alj|2002j ). In a 
smooth medium the dust temperature is a continuous and 
monotonic function of the distance from the central source. 
In a clumpy medium, instead, the non-illuminated side of 
the clouds, particularly those with high optical depth, will 
be in general much colder with respect to their illuminated 
side. The net result is that a whole range of temperatures 
can coexist within the torus at a given dis tance from the 
primary source (see Schartmann et al] |2008l . The tempera- 
ture of the dust within each individual cloud spans a wide 
range and therefore each region of the torus emits strongly 
at all wavelengths, from mid-IR to submm. The rms of the 
Wir distribution of smooth models is twice as large as that 
of clumpy models and smooth models can also produce large 
Wir. In fact, the largest ones in this study are indeed pro- 
duced by the smooth models in the matched grids. 

The peak of the IR SEDs, X pea k, measured on the con- 
tinuum alone, excluding the silicate feature is shown in the 
middle panel of Fig. [5] Both sets of models peak at wave- 
lengths typically between 10 and 30 um, a fraction of smooth 
models, however, have their peak at very short wavelengths. 
For both views, the majority of clumpy models culminate at 
~30 um, a behaviour related but not confined to the wider, 
on average, clumpy SEDs. 

The distribution of the near-IR spectral index, uir, 
varies a lot depending on the dust morphology (Fig. O lower 
panel). For type 1 views in particular, the range of overlap 
is very small, with clumpy models producing, on average, 
steeper spectra. The mean (median) of the distributions for 
the smooth and clumpy models are 0.63 (0.58) and 1.63 
(1.74), respectively. The same occurs for type 2 views but to 
a lesser extent, with a mean (median) of 2.98 (2.1) and 3.32 
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(3.18) for smooth and clumpy dust, respectively. Three com- 
ponents contribute to these differences: i) the lac k of a very 
hot c omponent in the clumpy models (see e.g. iDeo et all 
|2011| ). ii) the primary source (different power laws and dif- 
ferent wavelength coverage for the two dust distributions), 
and, to a much lesser extent, iii) the scattering in the clumpy 
medium, which strongly depends on the distribution of the 
clouds. We calculated the average fractional contribution of 
the primary source to the flux at 3.6 and 4.5 um to be of 
1.5% and 0.7%, respectively, for clumpy models, while it is 
constant and ~44% in smooth models in both bands. This 
would account for about 30% of the difference in the values 
of aiR between the two dust configurations and, therefore, 
the remaining difference must be due to point i. 

Some recent works indicate that the mid-IR emission, in 
particular the monochromatic luminosity at 12.3 um, L12.3, 
can be used as a diagnostic t o distinguish betwe en a smooth 
and a clumpy configuration. iHorst et alj (|2006T l found that 
in a sample of eight Seyfert galaxies (a Seyfert 1, an 1.2 
and six 1.5 or later), L12.3 is tightly correlated with their 
X-ray L2-iokeV luminosity, regardless of their type. This 
was interpreted as an evidence for the dust being optically 
thin at 12.3 um, a characteristic which was reported to be 
typical o f clumpy models, but incompatible with the smooth 
model of lPier fc Krolikl (l992^ . 

We have checked the sets of clumpy and smooth models 
against this prediction, by computing the ratio of L12.3 in 
type 2 over type 1 views, shown in Fig. [6] With the exception 
of 7% of clumpy models that are very close to be optically 



thin, i.e. with L[ y 2 p p jh\ 



1, both sets of models always 



lie well below this value. Overall, although the shape of the 
two distributions is different, we find no real evidence to the 
above claim. 

We have shown that even models with matched param- 
eters can produce very different features (see Fig. [3|. But 
could random parameter combinations result in very similar 
model SEDs? In order to answer this question, we introduce 
a measure of the dissimilarity of any two SEDs based on the 
properties discussed in the previous sections as follows: 
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(7) 



where sm and cl denote smooth and clumpy models, respec- 
tively, C = -F24/-F4.6 and F4.6 and F24, are the fluxes at 4.6 
and 24 /jm respectively. C has been introduced to take into 
account the relative fluxes of the (normalised) models in ad- 
dition to the IR spectral index, qjjj. Comparing the value 
of A for each pair of matched models, we find that none of 
the pairs have a value lower than 0.1. Furthermore, we only 
find 6 combinations for type 1 views and 21 combinations 
for type 2 views with A ^ 0.1, with very different model pa- 
rameters. We therefore conclude that there are virtually no 
random pairs of smooth and clumpy models in the matched 
grids that can produce very similar SEDs. 



4 DISCUSSION 

Within the paradigm of the AGN unified scheme, the dif- 
ferent properties of the various types of AGN are attributed 
to the interception of the line of sight by an axisymmetric 
molecular dust distribution encircling the nucleus, likely in 
the form of a torus. Direct and indirect observational ev- 
idence, such as the dichotomy of the AGN population in 
obscured and unobscured objects, the ionization cones for 
the NLR, the MIR bump present in the SEDs of all known 
AGN attributed to emission by dust heated by the primary 
source, but also IR interferometry that directly resolved the 
inner parsec of the nucleus of the pro t otype Seyfert 2 g alaxy 
NGC 1068 l|Wittkowski et all 12004 Ijaffe et al.1 l2004h , all 
support the existence of a toroidal absorbing structure. As 
there is evidence favouring both a clumpy d i stribution (e.g 
Risaliti et al. 1 120021 : IWittkowski et al.l 120041: I Tristram et al 



20071 ) and a smooth distribution (e.g. Ilbar fc Liral 1200 



both smooth and clumpy models are still equally widely 
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Figure 5. From top to bottom: W/ij, A pea j. and ck/a, for type 1 
(left column) and type 2 (right column) views, for smooth (con- 
tinuous lines) and clumpy (dashed lines) models. 



used to explain the observed SEDs of AGN, counting both 
successes and shortcomings. 

In this paper we compared two sets of models widely 
used in the literature, each representative of one of the two 
dust distri butions, namely an updated grid of the smooth 
models by Fritz et al.l {2006) and the clumpy models by 
iNenkova et al.1 (|2008a ) . In order to compare as similar mod- 
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Figure 6. The distribution of the ratio of £12.3 in type 2 over 
type 1 views, for smooth and clumpy models. 



els as possible and despite the intrinsic differences between 
the two configurations, we matched the two sets of param- 
eters building thus restricted grids, and compared only the 
models within these grids. This approach has the limitation 
of not exploring the full parameter space for either of the 
two models. Furthermore, and to avoid complications aris- 
ing from the probabilistic nature of viewing an AGN as a 
type 2 object through a clumpy medium, we consider two 
extreme inclinations i.e. edge on and face-on, for smooth 
models, and take as equivalent type 1 and 2 clumpy models 
those with a probability to see directly the central engine 
greater and lower than 0.5, respectively. Our findings can 
be summarised as follows: 

- Even after matching the model parameters there is 
not a one-to-one correspondence. For each smooth model 
of the restricted grid there can be several clumpy models 
with matched Y , q, and a, but different combinations of iVo 
and rv that finally produce a given value of T9.7. Addition- 
ally, and since we only consider the two extreme inclinations 
(face-on and edge-on), the torus opening angle can not be 
matched and for each clumpy model we consider all smooth 
counterparts, irrespective of the value of O. 

- The distribution of the various features of the IR model 
SEDs differs when smooth and clumpy dust configurations 
are considered. The very different behaviour of the silicate 
features is due more to the different chemical compositions 
assumed by F06 and N08 and less to the actual dust mor- 
phology. We confirm the occurrence of broader, on average, 
IR SEDs in clumpy configurations, with a larger fraction of 
clumpy models peaking at long wavelengths (~30 um). The 
infrared spectral index, especially for type 1 views, is the 
quantity that changes the most between the two dust con- 
figurations, owing to the differences in the primary source 
assumed in the two models as well as the lack of the hotter 
dust component in the treatment of the clumpy medium. 

- Our study showed no large differences in the behaviour 
of £i2 P 3 2 l^viz 1 1 as an indicator of how optically thin a 
medium is, between the two model configurations, with 5% 
of clumpy models showing a value very close to 1 (i.e. op- 
tically thin configuration), a similar amount of clumpy and 
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smooth models (~ 20%) having L*2 P 3 2 l^vz.i 1 ^ 0-9 an d the 
rest of both dust configurations lying below this value. 

- Models with matched parameters within the restricted 
model girds do not produce similar SEDs (similar either by 
eye or based on the value of A introduced in 33. 2[) . Addi- 
tionally, only a very limited number of random parameter 
combinations can result in seemingly identical SEDs, though 
the dust configuration differs. 

From the above we conclude that, even though the two 
dust models produce distinct SEDs, most of the differences 
arise from the model assumptions (e.g. primary source, dust 
chemical composition) and not from the dust morphology 
(smooth or clumpy). To summarise, the properties of dust 
in AGN as measured by matching observations (be it broad 
band IR photometry or 1R spectra) with models will strongly 
depend on the choice of the dust distribution. The possibil- 
ity to discriminate between a smooth and a clumpy medium 
based on the various SED features may exist, but ambigu- 
ities are more common than not. Independent estimates of 
physical parameters, such as the optical depth, the size of 
the torus or the mass of the gas are needed in order to further 
constrain the models. X-ray observations could, for instance, 
provide an upper limit of the optical depth, integrated along 
the line of sight, high resolution HI maps of known nearby 
AGN could put constraints on the gas content within the 
circumnuclear region. 

Eventually, ALMA will permit to indirectly determine 
the morphology of the obscuring material by allowing the 
comparison betwe en the mid-IR and sub -mm emission of the 
structure (see e.g. iMaiolino et al.| [2008) or even to directly 
resolve the obscuring torus, making use of its very high angu- 
lar resolution (sub-pc scale at the distances of nearby AGN 
at high frequencies). 



ACKNOWLEDGMENTS 

This work makes use of the lNenkova et al l (|2008al lbh models: 
http://www.pa.uky.edu/clumpy/! A.F. would like to thank 
Robert Nikutta and Moshe Elitzur for providing detailed 
explanation on their models and on the calculation of the 
output parameters, as well as Hagai Netzer for very useful 
comments. We also thank R. Nikutta for his details and 
insightful referee report, that improved a lot the content 
and quality of our work. This work makes use of TOPCAT, 
developed by M. Taylor. 



REFERENCES 

Agol E., Gogarten S. M., Gorjian V., & Kimball A., 2009, 

ApJ, 697, 1010 
Antonucci R., 1993, ARA&A, 31, 473 

Baes M., Verstappen J., De Looze I., Fritz J., Saftly W., 
Vidal-Perez E., Stalevski M., Valcke S., 2011, ApJS, 196, 
22 

Barvainis R., 1987, ApJ, 320, 537 
Berta S., et al., 2007, A&A, 467, 565 

Deo R. P., Richards G. T., Nikutta R., Elitzur M., Gal- 
lagher S. C, Ivezic Z., Hines D., 2011, ApJ, 729, 108 
Draine B. T. 2003, ApJ, 598, 1017 

Dullemond CP. & van Bemmel I.M., 2005., A&A, 436, 47 



Efstathiou A., & Rowan-Robinson M., 1995, MNRAS, 273, 
649 

Esquej P. et al., 2012, MNRAS, in press 
Fritz J., Franceschini A., Hatziminaogiou E., 2006, MN- 
RAS, 366, 767 
Granato G.L., Danese L., 1994, MNRAS, 268, 235 
Hatziminaogiou E. et al., 2010, A&A, 518, L33 
Hatziminaogiou E., Fritz J., & Jarrett T. H. 2009, MNRAS, 
399, 1206 

Hatziminaogiou E. et al., 2008, MNRAS, 386, 1252 
Hao L., et al., 2005, ApJL, 625, L75 

Hernan-Caballero A. & Hatziminaogiou E., 2011, MNRAS, 
414, 500 

Horst H., Smette A., Gandhi P., & Duschl W. J. 2006, 

A&A, 457, L17 
Ibar E., & Lira P. 2007, A&A, 466, 531 
Jaffe W., et al., 2004, Nature, 429, 47 
Krolik J.H & Begelman M.C., 1988, ApJ, 329, 702 
Laor A., & Draine B. T., 1993, ApJ, 402, 441 
Maiolino, R. 2008, New Astronomy Reviews, 52, 339 
Mason R. E., Levenson N. A., Shi Y., Packham C, Gorjian 

V., Cleary K., Rhee J., Werner M., 2009, ApJL, 693, L136 
Mathis J.S., Rumpl W., & Nordsieck K.H., 1977, ApJ, 217, 

425 

Mor R., Netzer H., & Elitzur M., 2009, ApJ, 705, 298 
Natale G., et al., 2010, ApJ, 725, 955 

NenkovaM., Sirocky M. M., Nikutta R., Ivezic Z., & Elitzur 

M. 2010, ApJ, 723, 1827 
NenkovaM., Sirocky M. M., Nikutta R., Ivezic Z., & Elitzur 

M. 2008b, ApJ, 685, 160 
Nenkova M., Sirocky M. M., Ivezic Z., & Elitzur M. 2008a, 

ApJ, 685, 147 
Nenkova M., Ivezic Z., & Elitzur M., 2002, ApJ, 570, 9 
Nikutta R., Elitzur M., & Lacy M. 2009, ApJ, 707, 1550 
Ossenkopf V., HenningT., & Mathis J. S. 1992, A&A, 261, 

567 

Pier E.A. & Krolik J.H., 1992, ApJ, 401, 99 

Polletta M., Weedman D., Honig S., Lonsdale C.J., Smith 

H., Houck J., 2008, ApJ, 675, 960 
Pozzi F. et al., 2010, A&A, 517, 11 
Ramos-Almeida C. et al., 2011, MNRAS, 417, L46 
Rees M. J., Silk J. I., Werner M. W., & Wickramasinghe 

N. C. 1969, Nature, 223, 788 
Risaliti G., Elvis M., Nicastro F., 2002, ApJ, 517, 234 
Rodighiero G. et al., 2007, A&A, 470, 21 
Rowan-Robinson M. 1995, MNRAS, 272, 737 
Schartmann M., Meisenheimer K., Camenzind M., Wolf S., 

Tristram K. R. W., & Henning T. 2008, A&A, 482, 67 
Schartmann M., Meisenheimer K., Camenzind M., Wolf S., 

& Henning T., 2005, A&A, 437, 861 
Shi Y. et al., 2006, ApJ, 653, 127 

Siebenmorgen R., Haas M., Kriigel E., Schulz B., 2005, 
A&A, 436, L5 

Sirocky M. M., Levenson N. A., Elitzur M., Spoon 

H. W. W., & Armus L., 2008, ApJ, 678, 729 
Stalevski M., Fritz J., Baes M., Nakos T., Popovic, L. C., 

2012, MNRAS, in press 
Sturm E., Hasinger G., Lehmann I., Mainieri V., Genzel 

R., Lehnert M. D., Lutz D., & Tacconi L. J., 2006, ApJ, 

642, 81 

Teplitz H. I., Armus L., Soifer B. T., et al., 2006, ApJL, 
638, LI 



© 0000 RAS, MNRAS 000, 



Smooth and Clumpy dust distribution in AGN 9 



Thompson G. D., Levenson N. A., Uddin S. A., & Sirocky 

M. M., 2009, ApJ, 697, 182 
Tristram, K. R. W., et al. 2007, A&A, 474, 837 
Urry, C. M. & Padovani P., 1995, PASP, 107, 803 
van Bemmel I. M., & Dullemond C. P. 2003, A&A, 404, 1 
Vignali C. et al. 2011, MNRAS, 416, 2068 
Vignali C. et al. 2009, MNRAS, 395, 2189 
Wittkowski M., Kervella P., Arsenault R., Paresce F., 

Beckert T., Weigelt W., 2004, A&A, 418, 39 



© 0000 RAS, MNRAS 000, 



